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The exchange interaction and correlations may produce a power-law decay instead of the usual 
exponential decrease of the wave function under potential barrier. The exchange- assisted tunneling 
vanishes in the classical limit, however, the dependence on the Planck constant h is different from 
that for a conventional single-particle tunneling. 
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I. INTRODUCTION 

In the Hartree-Fock equations the exchange interaction 
is described by the non-local (integration) operator K (r), 
and the well-known theorems proven for the Schrocdinger 
equation with a local potential U (r) are violated if we add 
the exchange term (or any other non-local operator). The 
Hartree-Fock equation for a fermion (e.g. bound atomic 
electron) at orbital \&&(r) has the following form: 



2m dr 2 



¥ 6 (r) + (£/(r) - E)9 b (r) = K(r) (1) 



(2) 



Here the summation runs over all occupied orbitals \& 9 (r) 
with the same spin projection as Vf'(r). The range of K(r) 
is determined by the range of the highest orbitals ^ , 9 (r) 
which stay outside the integral in Eq. |(2J). If ^q(r) be- 
longs to continuum (like an electron state in a conduct- 
ing band in a crystal), the range is infinite. In this case 
the range of the bound electron solution of Eq. {I) (e.g. 
atomic Is orbital) is also infinite since ^(r) can not van- 
ish if the right-hand-side of Eq. (Q]) is not zero. 

According to (HQ (see also HHH) the exchange can 
produce a power-law decay instead of the usual exponen- 
tial decrease in a classically forbidden region. For inner 
orbitals inside molecules decay is r~ 2 [lj, for macroscopic 
systems cos (kfr)r~ v , where kf is the Fermi momentum 
and v = 3 for ID, v =3.5 for 2D and v =4 for 3D crystal 
P, [H. Correlation corrections within the perturbation 
theory approach do not change these conclusions 0, [|[ . 
Slow decay increases the spin-spin interaction between 
localized spins in solids and the under-barrier tunnel- 
ing amplitudes. According to Ref. @ the exchange in- 
teraction may increase probability of ionization of inner 
atomic electrons by an external electric field (e.g. laser 
field) by many orders of magnitude (a different enhance- 
ment mechanism, an "atomic antenna", was suggested 
by Kuchiev [?J and rediscovered by Corkum Q). 

In this paper we want to consider the problem in 
the classical limit ft — ► 0. In the case of a single- 
particle problem a wave function under the potential 



barrier is exponentially small, ^ ~ exp (— J \p\dr/H), 
i.e. the probability of the tunneling tends to the clas- 
sical result (zero) very fast, it contains exp (—l/K). Here 
\p\ = \j2m(U (r) — E) is the semiclassical under-barrier 
momentum. At first glance, one may say that the ex- 
change interaction does not exist in the classical physics, 
therefore, the probability of the exchange-assisted tun- 
neling should also rapidly vanish. Moreover, any effects 
of identity of particles (Fermi and Bose statistics) should 
vanish when the wavelength A — h/p tends to zero. How- 
ever, the answer may be not as simple as it looks. For 
example, let us consider a three-body problem. 
A homogeneous electric field is applied to a hydrogen 
atom where the electron energy is E b . If the field is not 
too strong, the probability of the tunneling to the state 
of the same energy E b in continuum outside the atom is 
practically zero. 

Then another electron with energy E c (above the bar- 
rier) hits the atomic electron, loses energy E c — E b and 
occupies the state with energy Eh (equal to the energy of 
the bound electron) in continuum outside the atom. 
The lost energy E c — Ef, is transfered to the bound elec- 
tron which gets enough energy (E c ) to pass the barrier. 
Both electrons are free and have the same energies E b 
and E c . 

This process imitates the "tunneling" of a particle with 
the energy (for non-identical particles one more colli- 
sion with the reverse energy exchange is needed to imitate 
the "tunneling" ) . This picture looks like a purely classi- 
cal one. Therefore, it is not obvious why the exchange- 
assisted tunneling must vanish for Ti — > 0. A related 
question: if this effect vanishes, is the law exponential 
(similar to the vanishing of the single-particle tunneling) , 
or does it vanish less rapidly? 



II. TWO- WELL POTENTIAL 

To answer these questions it is instructive to consider 
a simple model with a minimal number of discrete states 
involved (to avoid interference of different contributions 
during summation and integration over spectrum which 
may change the law of Ti dependence; this happens for 
the dependence on distance r Q). 

We start from a model of resonance tunneling from 
one potential well to another potential well. The case 
of symmetric double-well potential has been solved e.g. 
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in the textbook Q. There are two levels corresponding 
to the symmetric (ground state) and antisymmetric wave 
functions. The tunneling produces the splitting of these 
levels, E± = E\ =F t\ where t\ ~ exp (— J \p\dr/h) is the 
tunneling amplitude, and the integral is taken between 
the classical turning points. 

If the first potential well ( "L" ) is slightly deeper than 
the second potential well ("R"), the ground state wave 
function may be presented as ip g = V'il + BtiipiR where 
Bti ~ t\/{E\L — Ei R ). Here we assume that the dis- 
tance to other levels is large, t-y <C (E±l ~ Eir) an d the 
probability of the particle in the ground state to be in 
the well R is exponentially small (proportional to the 
squared tunneling amplitude, B 2 X ~ t\/{E\L — Eir) 2 ). 

Now we add a second particle (identical fermion or bo- 
son) to a higher state 2 which has energy close to the 
top of the barrier. We can present its wave function as 
ip2 = ^2^2L + S2"02fl where the coefficient Bi is not nec- 
essarily small. In this case the probability of the particle 
in the ground orbital to be in the potential well R is no 
longer proportional to the exponentially small parameter 
t\. Indeed, the following two-step process takes place. 

Step 1: the second particle tunnels from the potential 
well L (orbital ip2L,) to the potential well R (orbital iI>2r)- 

Step 2: two-body process 2R, 1L — > 1R, 2L due to a 
nondiagonal Coulomb exchange interaction which trans- 
fers the first particle from orbital 1L to the orbital 2L 
and the second particle from 2R to 1R. 

As a result of these two steps, we have no change in the 
occupation of the state 2 and transfer of a particle from 
the ground state 1L to IR. This gives the amplitude for 
the ground state particle to be in the well "R" : 



B, 



Gl 



G(2,lL;li?,2) 
E\l — Em 



where 



G(2,1L;1R, 2) = 



(3) 



(4) 



i,2 {v)Hil (r) j- 6 — x iPirW ^2 (v')dv'dv 
r — r 



is the Coulomb exchange integral. Note that the poten- 
tial wells here may have one, two or three dimensions. 

This result may also be derived from the Hartree-Fock 
equations (|II2[) for the orbital Wf, — ipi — ipiL + Stpi, 
by projecting it to the orbital ipi R . Note that the con- 
tribution of the direct term in the Coulomb interaction 
between the particles 1 and 2 is included into the mean 
field potential U(r). 

Equation Q gives us dependence of the amplitude Bqi 
on the distance I between the wells L and R. If the dis- 
tance I ^> r\ where r% is the size of the orbital 1L, we can 
expand l/\r — r'\ near |r — r'| = /. Integral with the first 
term l/l of this expansion vanishes due to the orthogo- 
nality of the wave functions ipi(r) and ^(r). Therefore, 
the expansion starts from l/l 2 . The correlation effects 
correspond to higher orders in the perturbation theory 
in the Coulomb interaction integrals G, so they decay 
with distance faster than l/l 2 . 



III. CLASSICAL LIMIT 

Now consider the classical limit h — * 0. The exchange 
term in the Hartree-Fock equations (|H2p seems to ex- 
tend to a large distance (under the barrier) if we con- 
sider the exchange interaction of an under-barrier par- 
ticle (^b = "01 ) with an above-barrier one (^> q = fa)- 
The second term (the first non-vanishing term) in the 
multipole expansion of K(r) gives us ip2(r)/r 2 decay 
in the under-barrier area in the right-hand side of Eq. 
(fT]). Contrary to the single-particle tunneling amplitude 
(ti ~ exp (— J \p\dr/K)) this dependence does not con- 
tain h. How the classical limit may be reached in this 
case? The T% dependence actually comes from the coef- 
ficient before ip2(f) / 'r 2 which is given by the integral in 
K(r). 

The answer here depends on what kind of states we 
consider. 

1. Two semi-classical states, one below the barrier, an- 
other above the barrier. In this case we have integral be- 
tween two rapidly oscillating functions in K{r) in Eq.([2]), 
- J(dr'/|r - r'|) cos (J Pl {r)\dr/h) cos (J p 2 (r')dr' /ti). 
The integrals of this type are considered in the textbook 
Q in the chapter about matrix elements between semi- 
classical states. Such integrals are exponentially small 
at h — ► (they contain exp (— 1/ti)), therefore, the coeffi- 
cient before ip2(r) in K(r) in Eq.® in the right-hand-side 
of Eq. (fTJ) (and the exchange- induced Sipi) vanishes. 

2. Wave function of the first particle, tpi, is essen- 
tially a quantum state, a ground state or a low en- 
ergy state. Here the result depends on the proper- 
ties of the potential U(r). If this potential is smooth 
near the minimum, we may use the oscillator wave 
functions to describe the low-energy states. The ex- 
change integral in this case contains product of the 
oscillator function and the rapidly oscillating semi- 
classical wave function, ~ exp (— £, 2 /2) exp (ipx/ti) — 
exp (-(£ - g ) 2 /2) ex p (-p 2 / '{2muoh)) , 

where £ = x^J muj/h. The exchange integral is exponen- 
tially small again (see the last exponent). 

3. A singular attractive potential ( a model of a di- 
atomic molecule). In this case the size of the lowest 
states so rapidly tends to zero at h — > that there are 
no oscillations inside. Consider, for example, a Coulomb 
ground state where the size is given by the Bohr radius 
as = h 2 /(me 2 ). At r < as the product (p/h)r < 1, 
i.e. there are no oscillations. Therefore, K(r) has no 
exponential suppression for Ti — > 0. However, the size of 
the ground state (~ as) tends to zero and the binding 
energy(~ e 2 /as) to infinity, so this case actually has no 
classical limit. 



IV. CONCLUSION 

We see that the exchange-assisted tunneling vanishes 
in the classical limit (it still may be exponentially en- 
hanced in comparison with the single-particle tunneling) . 
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The difference with the classical "ionization by collision" 
case is that in the latter we do not fix final energy of 
each particle (provided the total final energy is equal 
to the initial energy, Eif + £?2/ = En + E^i), the en- 
ergy exchange depends on the initial conditions. In the 
exchange-assisted tunneling the spectrum of particle en- 
ergies in the final state is exactly the same as in the initial 



state (Eif — E2i, E%f — En), therefore, we, in fact, have 
imposed an additional constraint. For classical processes 
this event has zero phase volume (zero probability). 
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